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ABSTRACT 

The observed velocities of the gas in barred galaxies are a combination of the azimuthally- 
averaged circular velocity and non-circular motions, primarily caused by gas streaming along 
the bar. These non-circular flows must be accounted for before the observed velocities can be 
used in mass modeling. In this work, we examine the performance of the tilted-ring method 
and the DiskFit algorithm for transforming velocity maps of barred spiral galaxies into ro¬ 
tation curves (RCs) using simulated data. We And that the tilted-ring method, which does 
not account for streaming motions, under/over-estimates the circular motions when the bar 
is parallel/perpendicular to the projected major axis. DiskFit, which does include streaming 
motions, is limited to orientations where the bar is not-aligned with either the major or minor 
axis of the image. Therefore, we propose a method of correcting RCs based on numerical 
simulations of galaxies. We correct the RC derived from the tilted-ring method based on a nu¬ 
merical simulation of a galaxy with similar properties and projections as the observed galaxy. 
Using observations of NGC 3319, which has a bar aligned with the major axis, as a test case, 
we show that the inferred mass models from the uncorrected and corrected RCs are signifi¬ 
cantly different. These results show the importance of correcting for the non-circular motions 
and demonstrate that new methods of accounting for these motions are necessary as current 
methods fail for specific bar alignments. 


Key words: Cosmology: dark matter; galaxies: kinematics and dynamics; galaxies: individ¬ 
ual: NGC 3319 


1 INTRODUCTION 

The rotational velocities obtained from gas observations are one 
of the most commonly used tools to study the mass distribu¬ 
tions of spiral galaxies. It is often assumed that the gaseous 
motions are purely circular due to their low velocity disper¬ 
sions. However, the presence of a galactic bar will induce non¬ 
circular motions due to the streaming along the bar. Since bars ac¬ 
count for two-thirds of nearby galaxies (eg. de Vaucouleurs et al. 
1991; Eskridge et al. 2000; Knapen et al. 2000; Whyte et al. 
2002; Marinova & Jogee 2007; Menendez-Delmestre et al. 2007; 
Barazzaetal. 2008; Aguerrietal. 2009; Marino va et al. 2009; 
Mendez-Abreu et al. 2010; Masters et al. 2011), it is important to 
properly account for these motions. 

The study of the mass distribution is done by comparing 
the observed RC with the expected contributions from the gas, 
stars and dark matter, assuming axisymmetry, and that the gas is 
in circular motions (e.g. de Blok et al. 2008). For the dark mat¬ 
ter component, an empirical or theoretical density profile is used. 
There are several dark matter density profiles suggested in the 
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literature. The two most commonly used are the NFW profile 
(Navarro, Frenk & White 1996) derived from cosmological simu¬ 
lations with a central cusp and the isothermal (ISO) model (see 
e.g: Carignan 1985; Carignan & Freeman 1985) with a central flat 
core (see also: Einasto 1969; Burkert 1995; Kravtsov et al. 1998), a 
more general form of the dark matter density profile is also given by 
Blais-Ouellette et al. (2001). Alternative model such as the Modi¬ 
fied Newtonian Dynamics (MOND) (Milgrom 1983a,b,c) do not 
have a dark matter component. Milgrom (1983b) stipulated that 
there is no need for a dark matter halo to explain the flatness of 
galaxies’ rotation curves (RCs) if the law of gravity is modified 
below a critical acceleration ao. A comparison between ISO dark 
matter models and MOND models, for a sample of dwarf and spi¬ 
ral galaxies, is given in Randriamampandry & Carignan (2014) but 
a larger sample is needed before a final conclusion can be made. 

In order to distinguish between both these and other mod¬ 
els, RCs that represent the gravitational potential are necessary. 
RCs are often derived from observed velocity maps of gas us¬ 
ing the tilted ring method (Rogstad et al. 1974). This method di¬ 
vides the velocity map into co-centered rings and infers a variety 
of kinematic parameters, including the circular velocities. It has 
been implemented in the GIPSY (Allen et al. 2011) task ROTCUR 
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(Begeman 1989). Unfortunately, ROTCUR does not eorrect for the 
non-circular motions. There are other methods that have been de¬ 
veloped to account for these motions, including the publicly avail¬ 
able DiskFit code (Spekkens & Sellwood 2007). DiskFit creates 
physical models from the observed rotation velocities and extract 
the non-circular motions using a minimization technique (for 
more details see : ^). DiskFit has proven to be able to detect weak 
bars (Kuzio de Naray et al. 2012). It was specifically designed to fit 
a non-axisymmetric flow pattern to two-dimensional velocity maps 
(Spekkens & Sellwood 2007). However, it is well known that the 
DiskFit bar flow algorithm fails when the bar is almost parallel to 
the major or minor axis because of the degeneracy of the velocity 
components (Sellwood & Sanchez 2010). 

Dicaire et al. (2008) used numerical models with pure rotation 
and compared the known circular velocities with the observed ro¬ 
tational velocities. They found that if a bar is perpendicular to the 
major axis of a galaxy, the measured rotational velocities and thus 
the mass is over estimated and the contrary if the bar is parallel to 
the major axis. 

We propose an alternate method to account for non-circular 
motions using N-body simulations. A grid of models with differ¬ 
ent bar orientations, strengths, and bulge-disk-halo compositions 
is generated. RCs from images of these models are made using 
ROTCUR, which are then compared to the RCs calculated from the 
gravitational potential and its derivatives. This comparison gives 
the corrections necessary for non-circular motions. To get the cor¬ 
rected RC for a real galaxy, the best corresponding model from 
the grid of simulations is selected and the corrections derived from 
the simulation are applied to the actual galaxy’s RC derived from 
ROTCUR. 

The paper is organised as follow: Section 2 describes the de¬ 
tails of both the numerical simulations and the ROTCUR and Disk- 
Fit algorithms. Section 3 compares the RCs derived from the sim¬ 
ulations using ROTCUR and DiskFit to the RC obtained from the 
gravitational potential. In Section 4 we describe our new method of 
correcting RCs and use NGC 3319 as a test case to show the im¬ 
portance of properly accounting for non-circular motions in mass 
models. Finally, in Section 5 we present our conclusions. 


2 DATA ANALYSIS 

2.1 N-body simulations 

We use a series of spiral galaxy simulations, with different Hub¬ 
ble types, Sa, Sb, Sd, fitted to observed rotation velocities (cf 
Chilingarian et al. 2010) and different radial distributions of dark 
matter, with either a core or a cusp in the center, and different 
dark-to-visible mass ratios (maximum disk or not). They give a 
range of rotation velocities with various shapes, rising more or less 
slowly with radius. We do not vary the total mass of the galaxies, 
since only the shapes of the rotation velocities are important, and 
masses can be varied by a change of units. It is acknowledged that 
there is a mass and size evolution along the Hubble sequence (see 
e.g. van der Wei et al. 2014), and that late-type spiral galaxies are 
lighter and smaller, but these can be taken into account afterwards 
through rescaling. For all galaxy models, the total mass inside 35 
kpc is 2.4 10^^ M©, with 75% in dark matter and 25% in baryons, 
for the standard model (based on the ACDM cosmology) and equal¬ 
ity between baryons and dark matter for the maximum disk model. 

^ http://WWW.physics.rutgers.edu/-spekkens/diskfit/ 


We have used two different simulation codes, in order to test 
that our assumptions are not biasing the results. In both kinds of 
simulations, the runs are fully self-consistent, with live dark haloes. 
The first code is based on a 3D particle-mesh and the Fast-Fourier 
Transform (FFT) technique to solve the Poisson equation, and the 
gas is represented by sticky particles (FFT-sticky particles). We also 
ran a set of simulations using the Tree-SPH code used in the Galmer 
project (Di Matteo et al. 2007). More details on these simulations 
including the hydrodynamics, star formation, and feedback can be 
found in Appendix A. The total number of particles is 240k in the 
FFT-sticky code and 480k in the Tree-SPH code. 

These produced different outputs for the gas particle distribu¬ 
tions, but the effects of the bar on the observed velocity maps were 
the same. Since the conclusions are the same, we will primarily ex¬ 
amine the results for the FFT-stickly particles simulations and the 
results for the Tree-SPH runs can be found in Appendix B. 

For the initial conditions of the galaxies of different 
Hubble types, we have followed the procedure described by 
Di Matteo et al. (2007). The halo and the optional bulge are mod¬ 
eled as a Plummer sphere, and the stellar and gaseous disks follow 
a Miyamoto-Nagai density profile. 

The mass and characteristic radii of all components in the 
seven initial conditions, are displayed in Table 1. The standard 
model (-St) has a disk dominated by dark matter, with a central core. 
A NFW model, with a cusp and a central concentration of c = 8 has 
been run (-nfw), then a maximum disk model (-dmx) with equal 
mass between baryons and dark matter within the optical disk. The 
initial rotational velocities, corresponding to the main variations of 
the parameters, are displayed in figure 1. 

The expected velocities (Vexpected) were calculated from the 
gravitational potential and its derivatives. More specifically, the po¬ 
tential derived from the disk is averaged out azimuthally, to obtain 
T>(r), and the circular velocity Vdrc is 


2.2 Snapshots 

The snapshots from the FFT-sticky particles simulation were ex¬ 
tracted from the data using a projection code specifying the incli¬ 
nation angle and the position angle and deriving the radial velocity 
of the gas particles at each point of a data cube. The data cube have 
three axes, the first and second axis are the right ascension and dec¬ 
lination and the third axis is the velocity in km/s. Each cubes have 
50 channels with a 12.5 km/s velocity resolution. The cube is then 
transformed into a fits file. The projection angle PA was varied be¬ 
tween 10 to 180 degrees in order to cover a wide range of intrinsic 
bar orientations for all the snapshots. 

We select the snapshots where the bar has reached a quasi¬ 
stationary state (most of the gas have been converted into stars), 
but is still strong enough, and there is enough gas remaining in 
the disk to sample the rotation velocities with gas only. Indeed, the 
simulations do not take into account gas accretion form the cos¬ 
mic filaments, and the gas density is progressively reduced by star 
formation. 

2.3 Measuring the rotational velocities 

2.3.1 Using the tilted-ring method with ROTCUR 

ROTCUR uses the tilted-ring method (Rogstad et al. 1974), which 
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Table 1. Initial conditions for the various components (bulge, disk, halo, and gas) of the different models used. All models are normalized to giant galaxies to 
explore the different shapes. Dwarfs can be obtained through rescaling (see text). 


Model 

Mf, 

10“ M0 

kpc 

Md 

10“ M0 

kpc 

Mh 

10“ M0 

kpc 

10“ M0 

r^ras 

kpc 

gSa 

4.6 

2. 

6.9 

4. 

11.5 

10. 

1.3 

5. 

gSb-st 

1.1 

1. 

4.6 

5. 

17.2 

12. 

0.9 

6. 

gSb-dmx 

2.2 

1. 

9.2 

5. 

11.4 

12. 

0.9 

6. 

gSb-nfw 

1.1 

1. 

4.6 

5. 

17.2 

12. 

0.9 

6. 

gSd-st 

0. 

- 

5.7 

6. 

17.2 

15. 

1.7 

7. 

gSd-dmx 

0. 

- 

11.4 

6. 

11.4 

15. 

1.7 

7. 

gSd-nfw 

0. 


5.7 

6. 

17.2 

15. 

1.7 

7. 



> 




Figure 1. Initial RC s of the various models, described in table 1. The left plot compares the various morphological types, and the right plot the various dark 
matter distributions. 


divides the observed velocity field into concentric rings and es¬ 
timates the circular velocities and other kinematical parameters 
such as position angle (PA) and inclination using a minimiza¬ 
tion technique. This method assumes that the gas is moving purely 
along circular orbits. Each ring is characterized by the rotation cen¬ 
ter Xc and yc, the systemic velocity Vgys, the expansion velocity 
Vexp (assumed here to be zero), the inclination angle and PA. These 
parameters can be fixed or let free to vary depending on the user. 
The line of sight velocity at any position in the ring is then given 
by: 

V(x,y) = Vsys + Vcsin(i)cos(6>), (2) 

where 0 is the position angle with respect to the receding major axis 
measured in the plane of the sky, Vgys the systemic velocity and i 
the inclination angle. 

0 is related to the actual PA by the following equations: 
cosid) = -(^-^c)sin(PA) + (y-ye)cos(PA) 


“ Xc)cos(PA) + (y - yc)cos(PA) 

■ Rcos(i) 

The data cubes extracted from the simulation snapshots were 
analyzed with the Groningen Image Processing System (GIPSY: 
http : / /www. astro . rug . nl/-gipsy/) software. The task 
RFITS was used to convert the fits file into a GIPSY format, then 
the fits headers were changed using FIXHED. The following pa¬ 
rameters were modified in the header; the axis and the pixel size 
which is necessary before extracting the velocity fields. The task 
MOMNT was used to produce the moment maps, and the rotation 
velocities were derived using ROTCUR. The expansion velocity is 
fixed to zero for all the models. We choose an inclination angle of 
60 degrees which is suitable for the kinematical studies. The sys¬ 
temic velocity ( Vgyg) and position angle (PA) are fixed using the 
following values: Vgyg = 0.0 km s“^ and PA = 0.0 degree during 
the fitting procedures. 

The process of deriving the RCs includes three steps; firstly, 
the PA and inclination are held fixed using the optical values and the 
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systemic velocity and position centre (Xpos and Ypos) are allowed 
to vary; secondly, the systemic velocity and the kinematic centre are 
held fixed and the PA and inclination are let free to vary; finally, all 
the above parameters are fixed and the average values are used to 
derive the RCs. 


2.3.2 Using the DiskFit bisymmetric model 

DiskFit (Spekkens & Sellwood 2007) is specifically designed to 
correct the non-circular fiows due to bars. It creates physical mod¬ 
els from the observed velocity fields and extracts sthe velocities and 
other parameters from the model velocity map. DiskFit also uses a 
bootstrap method to estimate the uncertainty on the parameters and 
uses all to data during the fitting procedures. 

DiskFit compares the observed velocity field with the follow¬ 
ing equation: 

Vmodei = Vsys + sin(i) [VtCos(6>) - V2,tCos(26>b)cos(6>) 

(5) 

—V2,rsin(2^b)sin(^) 

where Vt is the circular velocity, V 2 ,t and V 2 ,r the amplitudes 
of the tangential and radial component of the noncircular motions 
for a bisymmetric flow model. This form of the equation is for 
bisymmetric fiows from bars. Moreover, DiskFit can also deal with 
warped galaxies. Regardless, if the gas is moving on purely circu¬ 
lar orbits, this equation reduces to the key equation of the tilted-ring 
method. 

DiskFit estimates the rotational velocities and the other pa¬ 
rameters using a minimization technique: 

^2 ^ A^ Vobs(x,y) -Ek=l^k,nVk ^2 

n=l 

where Vobs(x,y) is the observed velocity at the position (x,y) on 
the sky, Un is the uncertainty, ujk,n is a weighting function which 
includes the trigonometric factors and also defines an interpolation 
scheme for the projected model (Sellwood & Sanchez 2010). Some 
of the advantages of DiskFit are listed in Sellwood & Sanchez 
(2010). DiskFit uses all the data in a single fit, which is necessary 
to detect mild distortion and also gives a better estimate of the un¬ 
certainty using a bootstrap method (Spekkens & Sellwood 2007). 


3 RESULTS AND DISCUSSIONS 
3.1 Results using rotcur 

The results for the FFT-sticky particles code for three different bar 
orientations are shown in figures 2, 3 and 4 for the gSa, gSb and 
gSd models respectively. The snapshots were taken at the follow¬ 
ing epoch T=700 Myrs for gSa, T = 700 Myrs for gSb and T = 
900 Myrs for gSd. The intensity maps of the particles, which show 
the orientations of the bar (MOMENTO maps) are shown in the 
left panels and the velocity fields with iso-velocity contours (AV 
= 50 km s“^) (MOMENT 1 maps) are shown in the middle. The 
axes of the moment maps are in arcsec. A distance of 16.5 Mpc is 
adopted to convert from arcsec to parsec scale for all the models. 
Comparison between the expected RCs and the measured RCs us¬ 
ing ROTCUR are shown in the right panels. The expected RCs are 
presented as continuous lines and the RCs derived using ROTCUR 
as black circles. 


Table 2. Maximum departure from the average velocity (i.e maximum de¬ 
viation) for the snapshots with the strongest bars. 


Model 

Epochs 

Myr 

maximum deviation 

km 

gSa-st 

800 

60.72 

gSb-st 

1000 

51.45 

gSb-dmx 

600 

53.45 

gSb-nfw 

900 

52.09 

gSd-st 

1000 

57.27 

gSd-dmx 

700 

68.48 

gSd-nfw 

900 

25.80 


ROTCUR over-estimates the rotation when the bar is perpen¬ 
dicular to the major axis except for the gSa model for which most of 
the mass comes from the bulge component. The top panels of figure 
2 shows that the measured rotation velocities using ROTCUR and 
the expected rotation velocities are in good agreement even if the 
bar is perpendicular with the major axis. This is probably because 
the higher velocity due to the gas streaming motion along the bar 
is compensated by the contribution from the bulge. The differences 
between the measured rotation velocities and those expected from 
the gravitational potential are the largest when the bar is aligned 
with the major axis. Figure 5 shows a comparison between the ex¬ 
pected rotation velocities from the total mass of the galaxy and the 
average rotation velocities (obtained using ROTCUR) for four snap¬ 
shots. The rotation velocities of 12 simulated galaxies with differ¬ 
ent bar orientations were extracted and averaged for each snapshot 
(RA. varies between 0 and 180 in steps of 15 degrees). The er- 
rorbars are the 1-cr standard deviations. In the top-left panel, the 
bar starts to form but its effect on the derived rotation velocities is 
small. The two bottom panels present larger errorbars, with a max¬ 
imum of 60 km s“^ in the inner regions. 

The average RCs derived from the snapshots with the 
strongest bars (see Table 2 for the epochs) for all the models are 
shown in figure 6. The top panel is the gSa model with an isother¬ 
mal dark matter halo. The average of the measured RCs seems to be 
lower compared to the expected rotation velocities calculated from 
the gravitational potential. The errorbars are larger in the inner parts 
of the modeled galaxies which implies that the differences between 
the rotation velocities are only due to the bar orientations. These 
errorbars represent the variation of the derived RCs compared to 
the mean RCs. The gSb model (FFT-sticky) is shown in the mid¬ 
dle panel. On the left is the model with the isothermal dark matter 
halo, in the middle the maximum disk model and on the right the 
NFW model. The gSd model is presented in the bottom panel; the 
standard model on the left, the maximum disk model in the middle 
and the NFW model on the right. 

The maximum departure from the average rotation velocities 
(maximum standard deviation) for each model, are given in table 2. 
The first column is the name of the model, the second column the 
epoch corresponding to the strongest bar and the last column the 
maximum deviation from the mean velocity. 

3.2 Bisymmetric model results using DiskFit 

An input parameter file is used to run DiskFit. This file includes 
the following: a fits file that contains the velocities which are used 
to derive the kinematics; the initial values for the position angle 
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Figure 2. Three different bar positions for the model gSa (FFT-Sticky particles, T= 700 Myrs) Top panel: the bar is perpendicular to the major axis, middle 
panel: intermediate position and bottom panel: parallel to the major axis. The first column is the momentO map, second column the moment 1 map superposed 
with the iso-velocity contours and the third column the comparison between the expected velocities shown as a continuous line and the measured rotation 
velocities derived using ROTCUR (black points). The axes of the moment maps are in pixels (1 pixel = 5 arcsec). 


of the disk and bar, the systemic velocity, the axis ratio and the 
kinematic center. These parameters can be fixed or let free to vary. 
For a complete list of parameters the reader is advised to look at 
the original paper of Spekkens & Sellwood (2007) and the DiskFit 
manual. 

The results for the gSb model (T= 700 Myrs) using the bisym- 
metric model is shown in figure 7 for three different bar orien¬ 
tations. Top panel shows the bar perpendicular to the major axis, 
middle panel shows the bar in an intermediate position and on the 
bottom panel the bar is parallel to the major axis. The first column is 
the momentO map, second column the moment 1 map superposed 
with the iso-velocity contours and the third column the compari¬ 
son between the expected velocities shown as a continuous lines 
and the DiskFit results. The mean circular velocity Vt is shown as 


black circles, the tangential component of non-circular motion V 2 ,t 
as blue diamonds, radial components of the non-circular motion 
V 2 ,r as red squares and the expected circular velocity yexpected 
as a continuous magenta lines. T>b is the bar position angle in de¬ 
grees obtained from DiskFit. DiskFit is able to reproduce the ex¬ 
pected velocities calculated from the gravitational potential and its 
derivatives when the bar is in intermediate positions but runs into 
difficulties when the position angle of the bar $6 is parallel or per¬ 
pendicular to the major axis. This shows that DiskFit produces very 
large values for the rotation velocities and the errorbars when the 
bar is perpendicular with the major axis because of the degeneracy 
in the model parameters (Sellwood & Sanchez 2010). 

Figure 8 compares the projected position angle of the bar 
projected and thosc obtained using DiskFit. DiskFit gives higher 
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Figure 3. Three different bar positions for the model gSb (FFT-sticky particles , T = 700 Myrs see figure 2 for details). 


values as compared to the expected position angle of the bar and 
run into difficulty when the bar is close to 90 degrees. 

This is the reason why a different approach is needed when the 
bar is aligned with the major or minor axis. Here we propose to use 
N-body simulations of barred galaxies to estimate the non-circular 
motions. 


4 TEST CASE : NGC 3319 

DiskFit and ROTCUR are both very successful methods for deriving 
RC s from velocity maps. However, Section 3 uses simulations to 
show that both methods fail at recovering the true RC when a bar 
is aligned with the major or minor axis of an image. An alternate 
approach may be to use simulations to correct the RCs for non¬ 
circular motions. This will avoid many of the issues with either 
ROTCUR or DiskFit. We illustrate this method using NGC 3319, 


which has its bar aligned with the major axis, as a test case. NGC 
3319 is a gas-rich barred spiral galaxy classified as SBcd(rs) by 
de Vaucouleurs et al. (1991) at a distance of 14.3 Mpc (Sakai et al. 
1999). We compare the inferred mass profile for the uncorrected 
and corrected RCs and show that the conclusions about the dark 
halo structure vary greatly. While our method still needs refine¬ 
ment, this result clearly demonstrates the importance of taking into 
account non-circular motions. 


4.1 Observations and data reduction 

The Hi observations were obtained using the C configuration 
of the Very Large Array (VLA) for a total observing time of 
5.75h on-source. Details of the observations are presented in 
Moore & Gottesman (1998). The raw data were retrieved from the 
archive, edited and calibrated using the standard NRAO Astronom- 
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Figure 4. Three different bar positions for the model gSd ( FFT-sticky particles, T = 900 Myrs see figure 2 for details). 


ical Image Processing System (AIPS) packages. Data cubes were 
then obtained and cleaned with the AIPS task IMAGR using ro- 
bustO weighting, which optimizes sensitivity and spatial resolution. 
1031+561 is the phase and amplitude calibrator and 3C286 the flux 
calibrator. The C array data provide the proper spatial resolution 
for mass model analysis. The data cubes have 64 channels with 
10.33 km s“^ velocity resolution and 16 arcsecs angular resolu¬ 
tion, which corresponds to a linear physical resolution of ~ 1 kpc 
at the adopted distance of 14.3 Mpc. The moment maps were ex¬ 
tracted from the data cubes using the AIPS task MOMNT. Figure 
9 shows four different maps for NGC 3319. The WISE W1 image 
is shown at the top left. This band is essentially free from dust ex¬ 
tinction and gives the best representation of the stellar content of 
the galaxy. The moment 0 map (Hi gas distribution) is at the top 
right, the moment 1 map (velocity held) at the bottom left and the 
moment2 map (velocity dispersion) in the bottom right panel. 

The tilted ring method, implemented in the GIPSY task 


ROTCUR, is used to derive the RC . The momentl velocity map 
was used as input for ROTCUR. The tilted ring analysis was done in 
three steps: 

• First, the position angle (PA.) and inclination angle (i) were 
flxed using the values of 37° and 58° found by Moore & Gottesman 
(1998) in order to estimate the systemic velocity and the kinemat- 
ical centre. Our measured rotation centre and systemic velocity of 
742 zb 2 km s“^ are consistent with Moore & Gottesman (1998) 
(see their table 3). 

• Second, the systemic velocity and the rotation centre were 
kept flxed using the values obtained in the first step and the posi¬ 
tion angle and inclination were allowed to vary to investigate their 
variation as a function of radius. 

• Finally in the final step, the kinematical parameters were all 
flxed and only the rotation velocities were derived. 

The results from the tilted ring analysis are presented in figure 10. 
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Figure 5. Standard gSb galaxy model: comparison between the expected RCs calculated from the gravitational potential as continuous lines and the average 
RCs for different bar orientations as hlled circles with the 1-cr error bars. The RCs of 12 simulated galaxies with different bar orientations were extracted and 
averaged for each snapshot (PA varies between 0 and 180 degrees with a step of 15 degrees). The epochs of the snapshots are shown on the graph. 



Figure 8. Comparison between the projected bar position angles and those 
obtained from DiskFit. 


The variation of the inclination as a function of radius is shown in 
the top panel, the position angle in the middle panel and the RCs 
in the bottom panel. The results from the receding half is presented 


as red squares, those from the approaching half as blues circles and 
those obtained using both halves as black diamonds. 


4.2 Model-based correction to the RC 

As discussed in section 3, the RC is under-estimated when the bar 
is aligned with the major axis because of the gas streaming along 
the bar. Therefore, the effect of those non-circular motions need to 
be estimated, and a correction needs to be applied to the RC before 
it can be used for the mass models. 

We introduce the following step to correct the RCs for the 
non-circular motion induced by the bar: 

• The correction is derived from the differences between the ex¬ 
pected RC of the model, which is calculated from the true gravita¬ 
tional potential of the modeled galaxy, and the observed RC mea¬ 
sured using the tilted ring method (ROTCUR) 

• The correction is scaled in terms of velocity using the Vmax 

ratio, i.e. and in terms of radius using the ratio be- 

tween the disk scale lengths, i.e. (:5li£^2££i) before being used to 

^d,model 

correct the RC of an actual galaxy. 

• The modeled galaxy (model) and the actual galaxy (object) 
should have similar bar properties, the distance also have to be 
taken into account. 


The correction is done using the following 


^corrected — ^uncorrected COTTCCtZOTl X 


^max,object 
^max,model 


(7) 
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Figure 6. This figure shows the comparison between measured and calculated RCs for all the models listed in table 1 for the 
snapshots with the strongest bar (the epochs are given in table 2). Lines and symbols are the same as in figure 5. 


where correction — Vobs,model ^expected? ^obs,model IS the 
observed velocities obtain from ROTCUR, Yexpected the expected 
circular velocites; ycorrected and yuncorrected are the corrected 
and uncorrected rotation velocities of NGC 3319 respectively. 

To apply the corrections, the data from the model were nor¬ 
malized by the optical scale length in radius and by the maximum 
velocity in velocities using the steps described above. The data is 
then extrapolated to obtain the correction at a given radius. The 
distance of NGC 3319 is used to convert from angular scale into 
physical scale. Comparison between the corrected and uncorrected 
RC is shown in figure 11 . 

Despite the fact that NGC 3319 is classified as a SBcd galaxy, 
the gSb model presented at the bottom panel of figure 3 is used to 
derive the corrections. This is because the gSd model does not have 
a bulge component (see table 1) and that the bulge-to-disk (B/D) 
ratio of the gSb model obtained from the ratio between the bulge 
mass and disk mass given in table 1 M^/M^ = 0.239 is very similar 


to the B/D = 0.24 of NGC 3319 estimated from the 3.6 micron 
surface brightness profile. 


In this illustrated test case, we do not pretend to accurately 
correct for the non-circular motions induced by the bar. In the con¬ 
tinuation of this work, we intend to match even more accurately 
the simulated galaxies to the actual galaxies. However, we think 
that the gSb simulation used is close enough to NGC 3319 to show 
clearly how not taking into account the effect of bars can lead to 
deriving very different parameters for the derived dark matter com¬ 
ponent. We still have to investigate if it is the best way to scale the 
model to the actual galaxies. However, this is surely close enough 
to serve the purpose of this study. 
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Figure 7. DiskFit results using the bisymmetric model (m=2) for the gSb model (T = 700 Myrs). Top panel: the bar is perpendicular to the major axis, middle 
panel: intermediate position and bottom panel: parallel to the major axis. The first column is the momentO map, second column the moment 1 map superposed 
with the iso-velocity contours and the third column the comparison between the expected velocities shown as a continuous line and the and the DiskFit results. 
The mean circular velocity Vt is shown as black circles, the tangential component of non-circular motion \ 2 ,t as diamonds, radial components of the non¬ 
circular motion as squares and the expected circular velocity y expected as a continuous lines, is the bar position angle in degrees obtained from DiskFit. 
The axes of the moment maps are in pixels (1 pixel = 5 arcsec). 


4.3 Mass model 

4.3.1 Dark Matter Halo component 

The distribution of DM is characterized by a density profile. 
The observationally motivated pseudo-isothermal halo (ISO) and 
the cosmologically motivated Navarro, Frenk and White (NFW) 
(Navarro, Frenk & White 1996) profiles are widely used. The mass 
model is done by comparing the observed RCs with the quadratic 
sum of the luminous and dark matter components using a chi- 
squared minimization technique. The mass models were carried out 
using the GIPSY tasks ROTMAS and ROTMOD. 


The RC is given by: 


V?ot — Vgas + V* + Vhalo (8) 


where Vgas is the gas contribution, 14 the stars contribution and 
Vhaio the contribution of the dark matter component. 
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Figure 9. This figure shows several maps of NGC 3319. The top left panel shows the WISE W1 band (3.4 microns) image, 
where the bar and spiral arms are clearly seen; the momentO map (Hi distribution) is shown in the top right panel, the moment 1 
map (Hi velocity field) in the bottom left panel and the moment2 map (Hi velocity dispersion) in the bottom right corner. 


The pseudo-isothermal (ISO) Dark Matter (DM) Halo Model 

The pseudo-isothermal (ISO) dark matter halo is a core-dominated 
type of halo. The ISO density profile is given by: 

= ( 9 ) 


where po is the central density and Rc the core radius of the halo. 
The velocity contribution of a ISO halo is given by: 


The NFW DM Halo Model 

The NFW DM halo profile was derived from cosmological simula¬ 
tions and it is commonly accepted in the ACDM framework. The 
NFW halo density profile is described by: 


Pnpw(R) = :R 7 r~^ 
Rs Rs 


Viso(R) = 


47rGpoR?(l 


Rc 


atan( —)) 


( 10 ) 


where pi ^ 3iTo / (SttG) is the critical density for closure of the 
universe and Rs is a scale radius. The velocity contribution corre¬ 
sponding to this halo is given by: 


We can describe the steepness of the inner slope of the mass 
density profile with a power law p ~ r“. In the case of the ISO 
halo, where the inner density is an almost constant density core, 
a = 0. 


Vnfw(R) = V200 


ln(l + cx) — cx/ (1 + cx) 
x(ln(l + c) - c/(l + c))] 


( 12 ) 
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Figure 10. Tilted ring model results obtained with the GIPSY task ROTCUR. 
The variation of the inclination as a function of radius is shown in the top 
panel and the position angle (PA) variation in the middle panel. The RC 
presented in the bottom panel was obtained with all the kinematical param¬ 
eters, including PA and the inclination, kept fixed to their average values. 
The results from the receding side are presented as red squares, those from 
the approaching side as blue circles and from both sides as black diamonds. 
The lines correspond to PA= 37° and i = 58°. 


where V 200 is the velocity at the virial radius R 200 , c = R 200 /Rs 
gives the concentration parameter of the halo and x is defined as 
R/Rs. The relation between V 200 and R 200 is given by: 
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Figure 11. Comparison between the model-corrected and uncorrected RCs. 


the momentO map into concentrated rings and estimates the density 
for each ring. The output from ELLINT is then used in ROTMOD 
to calculate the gas contribution in the mass model, assuming an 
infinitely thin disk. The profile was multiplied by 4/3 to take into 
account the primordial Helium contribution. 


Stellar Contribution 

The SPITZER 3.6 microns surface brightness profile is used to cal¬ 
culate the stellar contribution. Near-Infrared probes most of the 
emission from the old stellar disk population. It is also less af¬ 
fected by dust and therefore represents the bulk of the stellar mass. 
The data were retrieved from the SPITZER archive. After the fore¬ 
ground stars were removed, the images were fitted with concentric 
ellipses using the ELLINT task in GIPSY. The results were then 
converted into a surface brightness profile in mag/arcsec^ using the 
method in Oh et al. (2008). 

The surface brightness profile is corrected for inclination and 
galactic extinction before being converted into mass density. The 
method used by Oh et al. (2008) was also adopted to convert the 
luminosity profile into a mass density profile. The density profile is 
given by: 

S[Mopc“^] = X (14) 

where T*'® is the stellar mass-to-light ratio in the 3.6 microns band, 
/is.6 the surface brightness profile and is a constant used for 
the conversion from mag/arcsec^ to Lo/pc^. ® is given by: 


V 200 — 


R 200 X Ho 

loo 


(13) 


where Ho is the Hubble constant taken as Ho 
72.0 km s~^ Mpc“^ (Hinshaw et al. 2009). 


4.3.2 Luminous matter Components 
Gas Contribution 

The contribution from the gas is estimated from the momentO maps. 
The task ELLINT is used to compute the Hi gas density profile using 
the kinematic parameters obtained with ROTCUR. ELLINT divides 


C®-® = M|® + 21.56 (15) 

where Mq® is the absolute magnitude of the Sun in the 3.6 microns 
band. Using the distance modulus formula and the distance to the 
Sun, Oh et al. (2008) found: 

M|® = m|® + 31.53 = 3.24 (16) 

The bulge-disk decomposition 

The surface brightness profile in figure 12 clearly shows that this 
galaxy has a bulge in the center, thus a bulge-disk decomposition is 
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Figure 12. Bulge-disk decomposition, the green line is the best ht for the 
disk, the blue for the bulge and the dashed magenta line the sum of the 
contribution from both components. 

needed to better constraint the contribution from the stellar compo¬ 
nent. 

The disk component was fitted with an exponential disk using 
the following equation 

AtW = ^*0 + 1.10857^, (17) 

where is the disk scale length and /xo the central surface 
brightness. The remaining profile is fitted with a Sersic profile by 
letting the Sersic index n, the effective radius Re and the effective 
surface brightness /Xe as free parameters. The bulge is therefore 
described by 

A((r) =/ie + + 1], (18) 

Ihe 

where bn = L9992n - 0.3271 for 0.5 < n < 10 (Capaccioli 1989) 
The best fit result shown in figure 12 gives R^ = 5.77 di 0.06 kpc, 
/xo = 18.66 zb 0.02 mag for the disk and n = 0.87 zb 0.02, Re = 2.04 
zb 0.05 kpc and /Xe = 21.74 db 0.31 mag for the bulge. 

Mass-to-light ratio 

The mass-to-light ratio (T*) for the IRACl band is normally de¬ 
rived using the J-K color (see Oh et al. 2008, de Blok et al. 2008). 
However, this method gives a higher mass-to-light ratio for the 
disk component than the bulge component which is unphysical, 
de Denus-Baillargeon et al. (2013) uses chemo-spectrophotometric 
galactic evolution (CSPE) models to constrain the mass-to-light ra¬ 
tio of nearby galaxies as function of radius. We use their results 
to estimate the mass-to-light of NGC 3319. NGC 3319 and NGC 
2403 have the same morphological type (SBcd) and similar maxi¬ 
mum velocity but NGC 3319 has a larger bulge compared to NGC 
2403. The values used for the mass model are given in table 3. 

4.4 Results 

The mass model results for the uncorrected and corrected RC s are 
presented in figure 13 for the ISO and NFW models. The observed 


Table 3. Results for the ISO and NFW dark matter halo models using the 
uncorrected and corrected RCs. 


Halo Model 

Params 

Uncorrected RC 

Corrected RC 

ISO 

PO 

5.55±0.99 

74.08±15.15 


Rc 

6.91±1.11 

1.34±0.16 


Td 

0.25 

0.25 


Tb 

0.40 

0.40 


Xr 

1.35 

0.75 

NFW 

C 

0.43±0.01 

5.18±0.73 


R 2 OO 

147.25±11.00 

56.46ib3.44 


Td 

0.25 

0.25 



0.40 

0.40 


Xr 

3.24 

1.44 


Notes: is the mass-to-light ratio for the stellar disk and T;, for the 

bulge. The central DM density po is given in units of 10“^ M© pc“^; Rc 
and R 200 are in kpc. 


RC is shown as filled black circles with error bars, the gas contri¬ 
bution as dashed green lines, the stellar disk contribution as dashed 
red lines, the stellar bulge contribution as long-dashed blue lines 
and the contribution from the dark matter halos as dashed magenta 
lines. The continuous lines are the best-fit models to the RC. Ta¬ 
ble 3 summarizes the results for the ISO and NFW models. The 
columns are described as follows, the first column is the type of 
DM halo used, the second column gives the parameters of the fits, 
the third and fourth columns are the results. For the ISO models, 
the dark matter halo is much more centrally concentrated when the 
corrected RC is used with a Rc ~ 1 kpc compared to ~ 5 kpc for 
the uncorrected RC and a much larger po (see table 3 for details). It 
is clear that the corrected and uncorrected RC give very different re¬ 
sults for the dark matter parameters which shows the importance of 
those corrections for non-circular motions. The same is true for the 
NFW models where the concentration index c is five times larger 
when the corrected RC is used and R 200 much smaller. Both the 
ISO and NFW models gives very good fits to the corrected RC 
with smaller Xr- This confirms that the parameters of the dark mat¬ 
ter halo are strongly affected by the non-circular motions induced 
by the bar, especially when the bar is aligned with the major axis. A 
large portion of the mass will be missed if the RC is not corrected 
and different conclusions could be drawn from the dark matter pa¬ 
rameters. The difference between the uncorrected and corrected RC 
is still significant even after the bulge component has been taken 
into account. This is a clear demonstration that corrections must be 
applied to the observed RC so that it really represents the dynami¬ 
cal mass of barred systems. 


5 SUMMARY 

We have presented the analysis of snapshot observations of N- 
body/hydrodynamic simulations of barred spiral galaxies with dif¬ 
ferent bar strengths and orientations to quantify the non-circular 
motions induced by a bar. . In the simulations, the bar forms natu¬ 
rally due to gravitational instabilities. Therefore, snapshots of sim¬ 
ulated galaxies spanning a wide range of bar strengths and bar ori¬ 
entations can be extracted and analyzed. Several models with dif- 
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Figure 13. Mass models using ISO DM halo model (top) and NFW DM halo model (bottom) for the uncorrected (left) and model-corrected (right) RCs. Lines 
and symbols are shown on the top left corner of the first panel. 


ferent morphological types and dark matter halo distributions are 
analysed. 

The results from the tilted ring method using ROTCUR show 
that the rotational velocities are largely under-estimated if the bar is 
parallel to the major axis and over-estimated if the bar is perpendic¬ 
ular except for an Sa galaxy which has an important central bulge 
for both simulations. This is consistent with Dicaire et al. (2008). 
This corresponds to the elongated orbits of the gas along the bar, 
which enhance the velocities at pericenter and reduce them at apoc- 
enter. Part of the gas and dust are moving alongside the bar which 
increases the measured circular velocity when the bar is parallel to 
the major axis and decreases it if the bar is aligned with the mi¬ 
nor axis. The effect is smaller when the bar is in an intermediate 
position. 

We also use DiskFit for our analysis. DiskFit is specifically de¬ 
signed to fit non-axisymmetric distortions from 2D velocity maps 
(Spekkens & Sellwood 2007). However, DiskFit works only when 


the bar is not aligned with the symmetry axes of the projeeted 
galaxy, and strong distortions are seen in the velocity field. 

Sellwood & Sanchez (2010) noted that the circular velocity is 
biased low when the bar is parallel to the major axis (apoeentre) 
since the gas follows an elliptical pattern and that the circular speed 
should be higher when the bar is aligned with the minor axis. 

The test case in section 4 shows that the distribution of mass 
in NGC 3319 is very affeeted by the non-circular motions induced 
by the bar. For example, the central dark matter core is much more 
compact when the corrected RC is used. The ISO and NFW results 
are very different in all cases for the corrected and uneorrected RCs. 

For future work, snapshot observations with different bar 
strengths and orientations will be studied to estimate the correc¬ 
tions that need to be applied to the observed rotation velocities of 
a sample of barred galaxies with a variety of bar strength and bar 
orientation. A sample of barred galaxies with different bar sizes 
and orientations will be seleeted and modeled. This will help us to 
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compare the results from the observations and the simulations for 
the individual galaxies.The main aim of this future work is to try 
to match as closely as possible the parameters of the model and the 
parameters of the galaxies, especially those of the bar. In particular, 
we will try to closely mateh the DM halo-to-luminous disk ratios. 
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APPENDIX 

A NUMERICAL DETAILS ABOUT THE SIMULATIONS 

The first simulation was done with a 3D Particle-Mesh code, based 
on FFT, with a useful grid of 128^. The algorithm of James (1977) 
is used to suppress the Fourier images, and gain a factor 8 in mem¬ 
ory. A total number of particles of 240 000 is distributed on the grid, 
and the density is computed with the cloud in cell algorithm (CIC). 
The gas is represented by sticky particles, colliding with inelastic 
collisions. At each encounter, the relative velocity of the clouds is 
multiplied by 0.65 in the radial direetion, and 1 in the tangential 
direetion, which keeps the angular momentum eonserved. 

The star formation is assumed to follow a Schmidt law, with 
exponent n=1.4, and a density threshold. It is ealibrated for a gas 
consumption time-scale of 5 Gyr. We do not assume instantaneous 
recycling, but take into aecount a continuous mass loss, across 
Gyrs (Jungwiert et al. 2001). The mass loss by stars is distributed 
through gas on neighboring particles, with a velocity dispersion of 
lOkm/s, to schematize the feedback energy. 

The RCs are fitted to the observedtypical curves for these 
Hubble types. For these giant spiral galaxies, the total mass inside 
35kpc is 2.4 10^^ M©, with 75% n DM and 25% baryons, for the 
standard model, and almost equality between baryons and DM for 
the disk-dominated (dmx) model. Aecording to the Hubble types, 
Sa, Sb, Sd, the bulge to disk ratio is varied. 

Similar simulations with TREE-SPH have also been run, 
with the same initial eonditions (c.f. Di Matteo et al. 2007). In 
this code, gravitational forces are ealculated using a hierarchical 
tree method (Barnes & Hut 1986) and gas evolution is followed 
by means of smoothed partiele hydrodynamies (SPH, Lucy 1977; 
Gingold & Monaghan 1982). Gravitational forces are ealculated 
using a toleranee parameter ^ = 0.7 and inelude terms up to the 
quadrupole order in the multipole expansion. A Plummer potential 
is used to soften gravitational forces, with same softening lengths 
for all partiele types, of 250pc. The eode and the various validation 
tests were described in Semelin & Combes (2002). In the present 
study, we use the isothermal gas phase (gas temperature of lO^K), 
adopted for the GalMer runs (Chilingarian et al. 2010). To capture 
shocks a conventional form of the artificial viscosity is used, with 
parameters a = 0.5 and /3 = 1.0 (Hernquist & Katz 1989). Be¬ 
cause of the short eooling time of disc gas, fluctuations in the gas 
temperature are quickly radiated away, so that simulations employ¬ 
ing an isothermal equation of state are justified. 

Star Eormation and continuous stellar mass loss, together with 


supernovae feedback are taken into aecount. We parametrized the 
star formation efficiency for a SPH particle as 
with the eonstant C chosen sueh that the unperturbed giant disc 
galaxies form stars with a gas depletion time-scale of 2 Gyr. This 
is equivalent to a Kennicutt-Sehmidt law of exponent 1.5. 

To avoid to form too small stellar particles, we adopt the 
method of hybrid particles, where progressively the gas particles 
eontain a stellar fraction for a certain time-scale, until the remain¬ 
ing gas fraetion falls below 5%. Then the particle is transferred to 
the stellar eomponent, and the gas is distributed over neighboring 
hybrid particles. 

The stellar feedback is implemented both with stellar mass 
loss, and through supernovae. It is assumed that a pereentage of 
1% of the young stellar population is becoming a supernova, i.e. all 
stars more massive than 8 M©, and the assumption of instantaneous 
recycling. The released mass also enriches the metallicity of the 
surrounding gas. 

The stellar mass-loss at each time step, is aceounted for by 
the method described by Jungwiert et al. (2001), and allows a more 
progressive metal enrichment of the interstellar medium. 


B RESULTS FROM THE SPH SIMULATION 

The dynamies of the gas in disk galaxies are very complex sys¬ 
tems and diffieult to reproduce in numerieal simulations. There¬ 
fore, to test the infiuenee of our assumptions on the gas treatment 
on our results, we have eonsidered the same simulations, with the 
same initial eonditions, but run with the TREE-SPH code, used 
in the GalMer data-base (cf Chilingarian et al. 2010). The main 
features of the SPH-tree code for this particular simulation and 
also our treatment of star formation and feedbaek, are described 
in the appendix A. The gas, which was considered essentially as 
an ensemble of dense clouds, dissipating their kinetic energy in 
collisions, is now considered as a fluid with pressure forees and 
shocks. The simulations have twice more particles, and the soften¬ 
ing length is twice smaller, so that the spatial resolution of the hy¬ 
drodynamic is larger. The velocity maps and other results from the 
GalMer project is publicly available and can be downloaded from 
their website (http : //galmer . obspm. fr/). Eigures B1 and 
B2 display the results with the TREE-SPH code for the gSa and gSb 
galaxies, corresponding to figure 2 and 3 with the EET-sticky parti¬ 
cles eode. The gas behaviour is different, as expeeted, but the global 
feature due to the bar, the perturbations of the RC, are similar. Eor 
the orientation where the bar is parallel to the line of sight, the dis- 
erepaney between the true circular velocity and the "observed" ro¬ 
tational veloeity is even larger, since the spatial resolution is higher. 
In the EET-sticky particles simulations, the high velocity gradient 
was a bit smoothed out. This is less visible, when the bar is per¬ 
pendicular to the line of sight, which is a configuration reducing 
the "observed" velocity gradient. Therefore, this shows that using 
different simulation eodes does not affect our results. 
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Figure Bl. Three different bar positions for the model gSa (SPH simulation, T= 200 Myrs). Top panel: the bar is perpendicular to the major axis, middle panel: 
intermediate position and bottom panel: parallel to the major axis. The first column is the momentO map, second column the moment 1 map superposed with 
the iso-velocity contours and the third column the comparison between the expected RCs shown as a continuous line and the measured rotational velocities 
derived using ROTCUR (black points). 
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Figure B2. Three different bar positions for the model gSb ( SPH simulation, T = 400 Myrs see figure B1 for details). 
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